Claims 

1 . A method for producing three-dimensional information of an 
object (4) in medical X-ray imaging, characterized in that 

the object is X-radiated from at least two different directions 
and the said X-radiation is detected to form projection data of the object (4) 

the object is modelled mathematically utilizing the projection 
data to solvej[s8] the imaging geometry and/or the motion of the object, where the 
solving concerns either some or all parts of the imaging geometry and/or the 
motion of the object. 

and said projection data and said mathematical modelling of 
the object are utilized in Bayesian inversion based on Bayes' formula 



pirn) 



to produce three-dimensional information of the object, the prior distribution ppr(9) 
representing the prior knowledge of the imaging geometry and/or the motion of the 
object (4), the prior distribution ppr(x) representing mathematical modelling of the 
object, x representing the object image vector, which comprises values of the X-ray 
attenuation coefficient inside the object, 0 representing the parameter vector of the 
imaging geometry and/or the motion of the object (4), m representing projection 
data, the likelihood distribution p(m|x,0 ) representing the X-radiation attenuation 
model between the object image vector x, geometry parameter vector 0, and 
projection data m, p(m) being a normalization constant and the posteriori 
distribution p(x,0|m) representing the three-dimensional information of the object 
(4) and the imaging geometry including the motion of the object. 
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2. A method according to claim 1, characterized in that the 
three-dimensional information of the object (4) is one or more two-dimensional 
images representing X-ray attenuation coefficient along slices through the object. 

3. A method according to claim 1, characterized in that the 
three-dimensional information of the object (4) is a three-dimensional voxel 
representation of the X-ray attenuation in the object. 

4. A method according to claim 1, characterized in that the 
measurement model is m = A9x+e, where matrix AG contains the weights 
describing how much each voxel contributes to the X-ray attenuation along the X- 
ray paths and the noise e is independent of object image vector x and the geometry 
parameter vector 0 leading to the likelihood distribution 

p( m \ x ,0) = p noise (m - A e x) . 

5. A method according to claim 1, characterized in that the said 
mathematical modelling comprises that X-radiation attenuates when passing the 
object (4), which means that every image voxel is nonnegative. 

6. A method according to claim 1, characterized in that 
mathematical modelling is expressed by the formula: 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the value 
U N (x) depends only on the values of voxels belonging to the neighborhood N, and 
a is a positive regularization parameter used to tune the width of the prior 
distribution. 
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7. A method according to claim 1 , characterized in that the 3D 
tomographic problem is divided into a stack of 2D tomographic problems and on 
every such 2D problem, the mathematical modelling is expressed by the formula: 

p pr (x) =exp^aY, U N (*)) 

N 

where the sum is taken over a collection of 2D neighbourhoods N and the value 
U N (x) depends only on the values of pixels belonging to the neighborhood N, and 
a is a positive regularization parameter used to tune the width of the prior 
distribution, and the 2D tomographic problems are related to each other by the 
formula 

P r3D(xG)) = exp(- y IZ | x(j)[k,q] - x(j-l)[k,q] |), 
where the sums are taken over all pixels (k=l,. . .,K, q=l,. . .,Q ) and y > 0 is another 
regularization parameter. 

8. A method according to claim 7, characterized in that the 
neighborhoods consist of two adjacent pixels and U calculates a power of the 
absolute value of the difference, leading to the formula 
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where s is a positive real number. 

9. A method according to claim 8, characterized in that s=l 
corresponding to total variation (TV) prior describing objects (4) consisting of 
different regions with well-defined boundaries. 

10. A method according to claim 1, characterized in that 
mathematical modelling is qualitative structural information of the target where the 
structural information is encoded in prior distributions that are concentrated around 
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object image vectors x that correspond to the physiological structures of the object 
(4). 

11. A method according to claim 1, characterized in that 
mathematical modelling consists of a list or probability distribution of possible 
attenuation coefficient values in the object (4). 

12. A method according to claim 1, characterized in that the X- 
ray imaging geometry, such as X-ray source position, has unknown error modelled 
in the distribution p(m|x,9). 

13. A method according to claim 1, characterized in that the X- 
radiation measurement noise is Poisson distributed. 

14. A method according to claim 1, characterized in that the 
medical X-ray imaging is dental radiography. 

15. A method according to claim 1, characterized in that the 
medical X-ray imaging is surgical C-arm imaging. 

16. A method according to claim 1, characterized in that the 
medical X-ray imaging is mammography. 

17. A method according to claim 1, characterized in that three- 
dimensional information of the object (4) is produced on the basis of the maximum 
a posteriori estimator (MAP) which is calculated by the equation: 

pOmap I m ) = max p(y I m ) 
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m representing projection data and y=(x,9) representing the vector composed of the 
object image vector and the geometry parameter vector and where the maximum on 
the right hand side of the equation is taken over all y. 



dimensional information of the object (4) is produced on the basis of the 
conditional mean estimator (CM), which is calculated by the equation: 



where m represents projection data and y=(x,0) represents the vector composed of 
the object image vector and the geometry parameter vector. 



dimensional information of an object (4) in a medical X-ray imaging, characterized 
in that the medical X-ray device (5) arrangement comprises: 

an X-ray source (2) for X-radiating the object from at least 
two different directions 

a detector (6) for detecting the X-radiation to form projection 

data of the object (4) 

means (15) for modelling the object (4) mathematically 
utilizing the projection data to solve the imaging geometry and/or the motion of the 
object, where the solving concerns either some or all parts of the imaging geometry 
and/or the motion of the object. 



means (15) for utilizing said projection data and said mathematical modelling of the 
object in Bayesian inversion based on Bayes' formula 



18. A method according to claim 1, characterized in that three 




19. A medical X-ray device (5) arrangement for producing three- 



and the medical X-ray device (5) arrangement comprises 



p(x,& | m) = 



Ppr(d)Ppr(x)P(™\x>6) 

p(m) 



25 



to produce three-dimensional information of the object, the prior distribution p pr (0) 
representing the prior knowledge of the imaging geometry and/or the motion of the 
object (4), the prior distribution p pr (x) representing mathematical modelling of the 
object, x representing the object image vector, which comprises values of the X-ray 
attenuation coefficient inside the object, 0 representing the parameter vector of the 
imaging geometry and/or the motion of the object (4), m representing projection 
data, the likelihood distribution p(m|x,G ) representing the X-radiation attenuation 
model between the object image vector x, geometry parameter vector 9, and 
projection data m, p(m) being a normalization constant and the posteriori 
distribution p(x,0|m) representing the three-dimensional information of the object 
(4) and the imaging geometry including the motion of the object. 

20. A medical x-ray device (5) arrangement according to claim 
19, characterized in that the three-dimensional information of the object (4) is one 
or more two-dimensional images representing X-ray attenuation coefficient along 
slices through the object. 

21. A medical x-ray device (5) arrangement according to claim 
19, characterized in that the three-dimensional information of the object (4) is a 
three-dimensional voxel representation of the X-ray attenuation in the object. 

22. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the measurement as 

m = A e x+e, 

where matrix A 0 contains the weights describing how much each voxel contributes 
to the X-ray attenuation along the X-ray paths and the noise e is independent of 
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object image vector x and the geometry parameter vector 0 leading to the likelihood 
distribution 

p(m \x,0) = p noise (m - A e x) . 

23. A medical X-ray device (5) arrangement according to claim 
19 characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically so that X-radiation attenuates 
when passing the object (4), which means that every image voxel is nonnegative. 

24. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by the formula: 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the value 
U N (x) depends only on the values of voxels belonging to the neighborhood N, and 
a is a positive regularization parameter used to tune the width of the prior 
distribution. 

25. A medical x-ray device (5) arrangement according to claim 
19, characterized in that the 3D tomographic problem is divided into a stack of 2D 
tomographic problems and on every such 2D problem, and the medical X-ray 
device arrangement comprises means (15) for modelling the object (4) 
mathematically by the formula: 
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p pr (x)=exp(raYU N (x)) 

N 

where the sum is taken over a collection of 2D neighbourhoods N and the value 
U N O) depends only on the values of pixels belonging to the neighborhood N, and 
a is a positive regularization parameter used to tune the width of the prior 
distribution, and the 2D tomographic problems are related to each other by the 
formula 

pr3D(xO» = exp(- y II | x(j)[k,q] - x(j-l)[k,q] |), 

where the sums are taken over all pixels (k=l,...,K, q=l,...,Q ) and y > 0 is another 
regularization parameter. 

26. A medical X-ray device (5) arrangement according to claim 

25, characterized in that the neighborhood systems consist of two neighboring 

pixels xj, xk or voxels xj, xk and Un ^ calculates a power of the 
/> r (^> )= ex P 

absolute value of the difference, leading to the formula where s is a positive real 
number and a is a regularization parameter used to tune the width of the prior 
distribution. 

27. A medical X-ray device (5) arrangement according to claim 

26, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by setting s=l corresponding to 
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total variation (TV) prior describing objects consisting of different regions with 
well-defined boundaries. 

28. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by assuming that mathematical 
modelling is qualitative structural information of the target where the structural 
information is encoded in prior distributions that are concentrated around image 
vectors x that correspond to the physiological structures of the target. 

29. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by assuming that mathematical 
modelling consists of a list of possible attenuation coefficient values in the object. 

30. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by assuming that the X-ray 
imaging geometry, such as X-ray source position, has unknown error modelled in 
the distribution p(m|x,0). 

31. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for modelling the object (4) mathematically by assuming that X-radiation 
measurement noise is Poisson distributed. 

32. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray imaging is dental radiography. 



29 



33. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray imaging is surgical C-arm imaging. 

34. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray imaging is mammography. 

35. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for producing three-dimensional information of the object (4) on the basis of 
the maximum a posteriori estimator (MAP), which is calculated by the equation: 

#Omap I m ) = maxp(y | m) 

m representing projection data and y=(x,0) representing the vector composed of the 
object image vector and the geometry parameter vector, and where the maximum 
on the right hand side of the equation is taken over all y. 

36. A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means 
(15) for producing three-dimensional information of the object (4) on the basis of 
the conditional mean estimator (CM), which is calculated by the equation 

^cm = jyp(y I m ) d y 

where m represents projection data and y=(x,0) represents the vector composed of 
the object image vector and the geometry parameter vector. 
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